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Abstract —Future machine to machine (M2M) communications 
need to support a massive number of devices communicating with 
each other with little or no human intervention. Random access 
techniques were originally proposed to enable M2M multiple 
access, but suffer from severe congestion and access delay in an 
M2M system with a large number of devices. In this paper, we 
propose a novel multiple access scheme for M2M communica¬ 
tions based on the capacity-approaching analog fountain code 
to efficiently minimize the access delay and satisfy the delay 
requirement for each device. This is achieved by allowing M2M 
devices to transmit at the same time on the same channel in 
an optimal probabilistic manner based on their individual delay 
requirements. Simulation results show that the proposed scheme 
achieves a near optimal rate performance and at the same time 
guarantees the delay requirements of the devices. We further 
propose a simple random access strategy and characterized 
the required overhead. Simulation resnlts show the proposed 
approach significantly outperforms the existing random access 
schemes currently used in long term evolution advanced (LTE- 
A) standard in terms of the access delay. 

Index Terms —Analog fonntain codes, machine-to-machine, 
message passing decoder, massive multiple access. 

I. Introduction 

M ACHINE-to-Machine communications have emerged 
as a promising technology to enable trillions of multi¬ 
role devices, namely machine-type communications (MTC) 
devices, to communicate with each other with little or no 
human intervention Qll. It has many potential applications, 
such as intelligent transportation systems (ITS), healthcare 
monitoring, retail, banking, smart grids, home automation and 
so on. It is expected that in the next a few years over 2 
billion MTC devices will become directly attached to cellular 
networks to provide M2M communications fsl]. Thus, there 
will be a massive number of MTC devices with no/low 
mobility [Ql in each cell, which is significantly more than 
the number of users in current cellular networks. Moreover, 
M2M traffic involves a large number of short-lived sessions, 
attempting to deliver a small amount of data (few hundred 
bits) to the base station, which is quite different from those 
in human-to-human (H2H) communications. Such differences 
motivate researchers around the globe to optimize the current 
cellular networks to effectively enable M2M communications 
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A. Motivation and Related Work 


In most existing wireless access networks, the first step in 
establishing an air interface connection, is to perform random 
access (RA) in a contention manner. In fact, short-lived 
sessions with small amount of data in M2M communications, 
makes it inefficient to establish dedicated resources for data 
transmission ||3l. Such a random access approach will not 
work effectively when the number of nodes is very large, 
due to frequent transmission collisions, leading to network 
congestion, unexpected delays, packet loss, high energy con¬ 
sumption, and radio resource wastage ||3l. Thus, one of critical 
challenges in M2M communications is to design an effective 
medium access scheme to support a large number of devices. 
In current RA approaches, when two or more devices select 
the same RA preamble in the first phase of the contention 
phase, a collision will occur and the respective devices will 
not be scheduled for data transmission. To reduce the access 
delay in M2M communications, several schemes have been 
proposed, such as dynamic allocation |^, slotted access, 
sroup-based O il llOT . pull-based, and access class barring 
techniques 111 lL|l2|] . Although these approaches can reduce the 
access collisions to a certain degree, most of them still suffer 
from very high access delays in highly dense networks. The 
main idea behind these schemes is to delay the retransmission 
of the access request for a random/fixed amount of time, thus 
increasing the access probability within a relatively short time. 
This is however inefficient in M2M communications due to 
small short burst transmissions of devices which mainly do not 
require the whole RB for their transmission. This means that 
for M2M communications with a very large number of devices 
there might not be enough RBs to be orthogonally allocated 
to the devices, which significantly increase the access delay 
even if the random access requests are delivered correctly. 


Additionally, different MTC devices have diverse service 
requirements and traffic patterns in M2M communications. 
Generally, we can divide the traffic types into four different 
categories. The first type is the alarm traffic, which is com¬ 
pletely random and its probability is very low; however, it has 
a very strict delay requirement. The second traffic type can be 
modeled by a random Poisson distribution with the parameters 
depending on the application 11131] . The regular traffic, such as 
smart metering applications, is the third traffic type, and the 
last type is the streaming, like video surveillance applications. 
Current proposals for enabling M2M communications, did not 
consider the priorities among devices and different quality 
of service (QoS) requirements. These approaches are mostly 
inefficient for M2M communications as they are generally 
designed for a fixed payload size and thus, cannot support 




2 


M2M applications with different service requirements. 

Recently, a systematic framework has been developed in 
EH to understand the fundamental limits of M2M commu¬ 
nications in terms of power efficiency and throughput. How¬ 
ever, they did not provide a systematic approach to develop 
an efficient communication protocol to approach these limits. 
Here, we consider a realistic model for M2M communications, 
which supports both regular and random traffics with different 
delay and service requirements. We develop a practical trans¬ 
mission scheme for M2M communications based on recently 
proposed analog fountain codes (AFCs) to enable massive 
number of devices communicating with a common base station 
(BS) while satisfying QoS requirements of all devices. We 
further show that the proposed scheme can closely approach 
the fundamental limits of M2M communications in terms 
of throughput and can satisfy the delay requirements of all 
devices at the same time. The main contributions of this paper 
are summarized next. 

B. Contributions 

1) Efficient Random Access Proposal: We consider the 
contention-based random access in the proposed probabilistic 
multiple access for M2M communications. In the contention 
process, each device transmits an RA preamble to the BS 
to establish a communication link with the BS. In existing 
RA schemes in M2M communications, the devices will be 
identified by the BS in the contention phase if each device 
has chosen an RA preamble which is not selected by other 
devices, and a resource block (RB) will be allocated to the 
device for the data transmission. For a system with a large 
number of devices, several devices may choose the same RA 
preamble and thus the BS cannot recognize the devices in the 
RA process. Moreover, even if the BS can identify all active 
devices, the number of RBs is not sufficient to support all of 
them. To overcome this problem, we propose a new contention 
mechanism. We group the RA preambles into several sub¬ 
sets according to delay requirements of MTC devices. Then, a 
sub-set of RA preambles is assigned for those devices which 
have the same delay requirement. In the contention period, 
each device selects a specific RA preamble based on its delay 
requirement from the associated subset of RA preambles. By 
detecting the RA preamble of the device, the BS then knows 
which subset this RA preamble belongs to and thus knows 
its delay requirement. The BS detects the number of devices 
which have selected the same RA preamble and broadcasts this 
information to the devices. The devices which have selected 
the same RA preamble then transmit in the same RB, and their 
identities will be transmitted along with their payload data and 
will be later recognized by the BS after the decoding. In this 
way, each available RB will be effectively used by a number 
of devices and RA collisions is handled in an efficient manner 
to simultaneously support collided devices. 

2 ) Efficient Maximum Throughput Transmission: Our sec¬ 
ond contribution in this work is to represent the multiple access 
process as a special kind of analog graph code. This enables us 
to draw on the recently proposed capacity approaching analog 
fountain codes iH to design the optimal access protocol 
achieving the sum-rate capacity of multiple access channel. 


Fountain codes have been used to facilitate the multiple 


access transmissions over erasure channels 11161 [17 


multiple 

HH, 


where collided packets are used in an effective way to increase 
the overall system throughput; however, they cannot be used 
for wireless channels which are affected by fading and noise. 
In lEoll . it has been shown that the capacity of a Gaussian 
multiple access channel can be achieved by using spatially 
coupled sparse graph multi-user modulation with an iterative 
interference cancelation scheme at the receiver. However, 
this approach performs decoding and interference cancelation 
separately and its complexity increases significantly as the 
number of access users increases. Thus, it will not be fea¬ 
sible for M2M applications with a massive number of users. 
To overcome this problem, we proposed a simple capacity- 
approaching analog fountain code (AFC) in jlSll with linear 
complexity in both the encoder and decoder. Thanks to the 
interesting linear property of AFCs, the summation of AFC 
coded signals of multiple simultaneous transmitted users still 
forms an AFC code 112111 . Thus the receiver can represent the 
multiple access transmission by a special kind of graph codes, 
and the BS can jointly decode them by using a standard belief 
propagation (BP) decoder. The proposed probabilistic multiple 
access AFC (MA-AFC) scheme enables the devices to transmit 
in a probabilistic manner in a way to form an equivalent AFC 
code at the BS, and ap proa ch the sum-rate capacity of the 
multiple access channel ll2l|] . 

3) QoS aware Random Transmission: Different devices in 
M2M communications have diverse traffic and service require¬ 
ments. Existing random access schemes do not provide QoS 
guarantees in the random access transmissions. This precludes 
use in many practical M2M systems, where timing constraints 
are critical. How to alleviate access congestions in massive 
access transmissions while providing delay or QoS guarantees 
remains a significant hurdle in M2M access network design. 
In the proposed approach, the degree distribution functions 
of AFCs and access probabilities are optimized such that the 
delay requirement for each device is satisfied. We formulate an 
optimization problem based on the density evolution analysis 
of the AFC code to find the optimum degree and access prob¬ 
ability for each device to satisfy their delay requirements. We 
further show that the proposed approach can simultaneously 
satisfy the delay requirement for all devices. 

Throughout the paper the following notations are used. We 
use a boldface small letter to denote a vector and the entry 
of vector v is denoted by rij. A boldface capital letter is used 
to denote a matrix and the element in the row of matrix 
G is denoted by g^ j. The transpose of matrix G is denoted by 
G^ A set is denoted by a curly capital letter and the cardinality 
of set S is represented by |5|. 

The rest of the paper is organized as follows. In Section 
II, we present the system model. The proposed probabilistic 
multiple access scheme for M2M communications is presented 
in Section III. Section IV presents the asymptotic analysis of 
the belief propagation decoding of the MA-AFC scheme based 
on the density evolution analysis. In Section V, an optimization 
problem is formulated to find the optimum code parameters. 
Simulation results are shown in Section VI, and finally Section 
VII concludes the paper. 
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Fig. 1. A resource block (RB) of duration Ts and bandwidth W. 

II. System Model 

In this section, we first introduce the network and channel 
model in M2M communications. Then we discuss is more 
details the packet transmission and delay model for MTC 
devices. 


A. Network Model and Random Access 

We consider an M2M last mile access system, where N 
MTC devices are uniformly distributed in a single cell of 
radius rg and communicate with a base station (BS) located at 
the origin. In this paper, we focus on the M2M applications 
where MTC devices are deployed at fixed locations, such as 
smart meters, sensors or cameras in buildings, roads, bridges, 
etc. 

It is assumed that the channel between each device and the 
BS is a slow time-varying block fading channel, for which the 
channel remains constant within one transmission block but 
varies slowly from one block to the other. We consider a time 
division duplex (TDD)-based wireless access system, where 
the channel gain of the uplink is assumed to be the same as 
that of the downlink Q. With this assumption, each device 
can estimate the uplink channel gain from the pilot signal sent 
periodically over the downlink channel by the BS. Using the 
pilot signal, MTC devices also synchronize their timing to that 
of the BS. The BS, however, does not have the knowledge 
of any channel state information (CSI). This assumption is 
particularly relevant in M2M communications, where due to a 
large number of devices, it would be impractical for the BS to 
obtain CSI to every MTC device 11241] . The channel between 
each device and the BS is modeled by path loss, large scale 
shadowing and small scale fading effects. Thus, the received 
power at the BS for the signal transmitted from is given 
by: 

Pr,r = ( 1 ) 

where j is the transmit power, 7 is the reference signal 
to noise ratio (SNR), defined as the average received SNR 
from a device transmitting at maximum power located at 
the cell edge, hi is the effective channel gain, defined as 


l^il = XiVi{n/ro)~ 


is the distance from to the BS, 


Xi is a log-normal random variable modeling shadowing gain 
with standard deviation a dB, r]i is small scale fading modeled 
by exponential random variable with parameter 1 , and a is the 
propagation loss factor. We further assume that the devices 
perform uplink power control such that the average received 
SNR from all devices at the BS are the same Ml311 and equals to 
7 o. Such an assumption has been considered for code division 
multiple access (CDMA) in IIMl . 

Time is divided into slots of duration ts secs and total 
available bandwidth is W Hz. Each time slot is referred to 


Notation 

Description 


Cell radius 

W 

Total bandwidth in Hz 

Ts 

The duration of a resource block (RB) 

T 

The ratio of the duration of RA process in an RB 

hi 

The effective channel gain of device Ui 

7 

The reference SNR 

7o 

Received SNR at the BS after power control 

A 

Average number of active devices in an RB 

Ai 

Packet arrival rate at the device 

N 

Total number of devices 

Ns 

Total Number of RA preambles 

k 

Payload size of each MTC device 


Vector of information symbols of Uj 


Generator matrix of the AFC code used in Uj 

dc 

AFC code degree 

Ws 

Weight set of AFC 

D 

AFC weight set size (AFC maximum degree) 

<^1 

Average weight power of an AFC 

Nd 

Number of delay groups 

ti 

The delay constraint of the 2 *^ delay group 

Pi 

Access probability of the device 

Lis.) 

Length of the random seed for devices with RA preamble s 

e 

The probability that BS fails to detect random seeds 


{ki, ki + 1, ... , ^2 } G Z for ki < k^ 


as a resource block (RB), where a fraction t, 0 < r < 1 , of the 
resource block is reserved for the random access procedure as 
shown in Fig. [T] Similar to lll3n . we assume that the number 
of active devices in each RB of duration ts is random and 
follows a Poisson process with rate A. We also consider a 
contention-based random access strategy in M2M communica¬ 
tions, where Ns different random access preambles, which are 
orthogonal frequency division multiplexing (OFDM) symbols 
in the FTE-Advanced standard ||3l, are selected in a random 
access attempt. The BS also assign an access probability to 
each device according to its delay requirement. The access 
probabilities will be later optimized in this paper to satisfy 
the delay requirement of all devices. The notation used in this 
paper is summarized in Table |I] for quick reference. 

It is important to note that in many M2M applications, M2M 
and human-to-human (H2H) communication coexists. H2H 
users have high priority to obtain a connection to transmit 
their data to the BS. Most of research work has considered 
that the contention-free random access is used for H2H users, 
that is the BS assigns a preamble and an RB to the H2H user 
in a timely manner lE^ . In this work we only focus on M2M 
devices, where opportunistically contend for RBs through a 
contention-based random access, and H2H users are assumed 
to have access to the BS through the allocated RBs. 


B. Delay Model for MTC devices 

In M2M communications, timing constraints for MTC de¬ 
vices typically range from 10 ms to several minutes due to the 
infrequent transmission features |[ll]. For instance, some M2M 
multimedia applications such as video surveillance systems 
with strict timing constraints have the delay requirement 
ranging from 10 ms to 40 ms. Such a diverse QoS requirements 
in M2M communications intensify the need for a practical 
channel assignment scheme to satisfy these requirements. 
Moreover, different MTC devices have diverse traffic patterns; 
i.e., some of them are regular and some have completely 
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random traffics. In Table [Hi we divide different M2M appli¬ 
cations based on their delay requirements and traffic patterns, 
according to the 3rd Generation Partnership Project (3GPP) 
M2M standard in ll]]. More specihcally, we divide the MTC 
devices into = 10 groups based on their applications, where 
the maximum affordable delay in group i is denoted by 

It is important to note that MTC devices with regular traffic 
have predictable activities and can be completely coordinated 
by the BS. This means that the BS will regularly assign an 
RB to the devices with regular activities and they can transmit 
their packets through the assigned RB. This can be easily ob¬ 
tained through several approaches in current LTE and LTE-A 
standards ll22ll . However, the devices with completely random 
activities should contend for RBs as allocating predetermined 
RBs to them is not efficient 111]. Therefore in this paper, 
we only consider MTC devices with random activities with 
irregular traffic patterns. This is the most challenging issue in 
M2M communications due to the ever increasing number of 
such devices. 

Each active MTC device U^, i e is assumed to have k 
information symbols to transmit within a maximum delay of 
tj, where l < j < N^. That is, these k information symbols 
must be delivered to the BS within tj seconds after 11^ starts 
transmitting this message, i.e, 11^ belongs to the delay 
group. We also assume that the device has the packet arrival 
rate A^, that is in each RB, packets are generated at the 
device, where each packet has to be delivered to the BS no 
later than tj seconds after its generation. Let us assume that 
l/Aj > tj. Therefore, if a previous packet has been delivered at 
the BS no later than tj seconds after its generation at the 
device, the new packet will be immediately sent by the device 
and there is no need to buffer it. On the other hand, if l/Aj < tj, 
we set the delay requirement of the device to l/Aj to avoid any 
further delay due to queuing. In other words, if the device has 
no memory to buffer the packets (i.e. each packet has to be 
successfully delivered at the BS before the generation of the 
next packet), the appropriate delay requirement for the device 
is min{l/Aj,tj}. It is also important to note that for devices 
with limited buffer size, the device can change the access 
probability according to the number of packets in its buffer. 
Moreover, a packet in M2M communications has a small size 
(few bytes); therefore, we can assume that all the packets in 
the buffer of a MTC device can be transmitted in one RB 
Q. Therefore In the rest of the paper, we only use the term 
delay requirement for each MTC device, and assume that the 
queuing delay has been taken into account when calculating 
the appropriate delay requirement of each device. 


III. The Proposed Rateless Probabilistic Multiple 
Access eor M2M Communications 

In this section, we develop an efficient multiple access 
scheme for M2M communications based on analog fountain 
codes. Eor this aim, we first briefly introduce AEC codes. 
Then, we present the rateless probabilistic multiple access for 
M2M communications, referred to as multiple access analog 
fountain coding (MA-AEC). 


TABLE II 

Several M2M applications with their delay constraints . 


Group 

Delay 

Constraint 

Example of 

application 

Traffic Type 

Message Size 

1 

10 ms 

Emergency 

Alarm 

Very Unlikely 

Very small 

2 

20 ms 

Intelligent 
transport 
system (ITS) 

Regular/lnegulai 

Medium 

3 

40 ms 

Video stream¬ 
ing 

Regular 

(Streaming) 

Laige 

4 

100 ms 

Control 

messages 

for devices 

o 

Regular 

Very small 

5 

500 ms 

eHealth 

Regular 

Medium 

6 

1 s 

Monitoring the 
Devices 

Regular 

Small 

7 

10 s 

eHealth 

Regular/lrregulai 

Medium 

8 

100 s 

eHealth 

Regular/Irregulai 

Medium 

9 

500 s 

Smart meters 

Regular 

Medium 

10 

>1000 s 

Smart meters 

Regular 

Medium 


A. Analog Fountain Codes 

Analog fountain codes (AEC) were originally proposed 
in lll5l l27n as an effective adaptive transmission scheme to 
approach the capacity of wireless channels. In AEC, the entire 
message of length k binary symbols is hrst modulated using 
binary phase shift keying (BPSK) modulation to obtain k 
modulated information symbols, bj e {-1,1}, where i e iJl- 
Then according to the predefined degree distribution function, 
dc randomly selected modulated information symbols are 
linearly combined with real weight coefficients to generate 
one coded symbol. Eor simplicity, we assume that the code 
degree, dc is fixed and weight coefficients are chosen from a 
finite weight set, Ws with D positive real members, as follows: 


Ws = {Wi e K+|i e zf} 


( 2 ) 


where R+ is the set of positive real numbers. Let us assume 
that b is a BPSK modulated vector of dimension fc x 1, and G 
is the generator matrix of dimension mx k, then AEC coded 
symbols, c, are generated as follows: 


c = Gb, 


( 3 ) 


where m is the number of coded symbols, only dc elements of 
each row of matrix G are nonzero, and each nonzero element 
of G is randomly chosen from the weight set Ws- To further 
enhance the performance, in Q we proposed to use a high- 
rate precoder to encode the original data before applying the 
AEC code. To fully utilize the constellation plane, i.e., both in 
phase and quadrature phases, two consecutive coded symbols 
compose one signal as We also use a standard 

belief propagation (BP) decoder to decode AEC codes. Eurther 
details of AEC encoding and decoding can be found in 1115 , 

il. 


B. Multiple Access AFC for M2M communications 

The proposed probabilistic multiple access AEC (MA-AEC) 
scheme contains two steps. In the first step, namely contention 
phase, the random access requests are sent by the devices to 
the BS and a connection is established between the devices 
and the BS. In the second phase, namely data transmission 
phase, the actual payload data along with the device identity 
(ID) is transmitted by the devices to the BS. 
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1) The contention phase (CP): To apply AFC codes to 
M2M communication systems and satisfy delay requirements 
of all devices, we need to slightly modify the random access 
procedure in LTE-A ||3l- The BS hrst partitions Ns RA 
preambles into < Ns subsets, according to delay constraints 
of MTC devices, and each group of RA preambles are assigned 
for one delay group. The BS then broadcasts this information 
to all devices. Let denote the set of RA preambles allocated 
to the delay group, then an MTC device with delay 
constraint will randomly select an RA from Si in the 
contention period. By detecting the RA preamble of the device, 
the BS will know which subset this RA preamble belongs to 
and thus know its delay requirement. 

In the hrst step of CP, an MTC device with delay constraint 
ti randomly selects an RA preamble from set Si and transmits 
it to the BS. Let s e Si denote the selected RA preamble, 
then, it is possible that more than one devices with the same 
delay constraint ti have chosen s. As we assume that all MTC 
devices perform power control, such that the received signals 
from different devices will have the same received power, 70 , 
the received power at the BS for the RA preamble s will be 
proportional to the number of devices which have selected RA 
preamble s. Let U{s) denote the set of MTC devices, which 
have selected the same RA preamble s. Since we assume 
that the RA preambles are orthogonal, the received signal 
corresponding to RA preamble s is given by: 

y = |W{s)|V^ s + z, (4) 


where z is the additive white Gaussian noise (AWGN) vector. 
The average received signal power at the BS can then be 
calculated as |W(s)|^ 7 o. Since the BS knows 70 , it can easily 
estimate |W(s)|. 

In the second step of CP, the BS transmits a random access 
response (RAR) message including the number of devices 
which have selected each RA preamble, the allocated RB 
for each detected RA preamble, and timing information and 
optimized access probabilities. In the third step of CP, each 
device hrst calculates the length of a random seed based on 
the received information from the BS, and then generates a 
random seed with the determined length and sends it to the BS. 
We assume that the random seeds are orthogonal sequences. 
Once the BS and each device have shared the same seed for 
the random encoding structure of ALC, they can construct 
the same bipartite graph to perform the ALC encoding and 
decoding. Pig. |2] shows the proposed random access procedure 
for M2M communications. 


Let L(s) denote the length of the random seed selected by 
devices which have selected the RA preamble s in the hrst step 
of CP. The BS fails to uniquely detect these random seede with 
the probability at most e, which is given by 10 : 


e 



2-£'(s) - 1 




( 5 ) 


This is followed from the fact that each random seed can be 
detected at the BS if it has been generated by only one device. 
Thus, the minimum length of the random seed L{s) required 
to have a non-detectable probability lower than a predehned 



Device BS 


Fig. 2. The proposed random access procedure for M2M communications. 


value e > 0 is as follows: 


L(S) 


log2 



1 — (1 — e) |W(s)|-l 


(6) 


where [.] is the ceil operator. 

In the proposed RA scheme, even when two or more devices 
select the same preamble, they will be allocated with the same 
RB, thus they are transmitting in the same RB. This means 
that the proposed approach can support many users with a 
negligible collision probability, which can be controlled by 
the length of the random seed. Unlike the conventional RA 
schemes, the devices in our approach can send their random 
access requests whenever they have data to transmit and do 
not need to wait for several RA attempt to get access to 
the network. In fact, once the device sends a random access 
preamble to the BS, it will be allocated with an RB, no matter 
whether that preamble has been selected by other devices or 
not. This will signihcantly decrease the access delay. 

2) The data transmission phase: After the contention 
phase, each active device and the BS have shared the same 
random seed. The devices which have selected the same RA 
preamble, are transmitting their ALC coded symbols in the 
same RB, in a probabilistic manner with the same access 
probability, but with different random seeds. Let pj denote 
the access probability assigned to Uj. It is important to note 
that devices which have selected the same RA preamble will 
be assigned with the same access probability. In the data 
transmission phase, in each time instant e, device U_, generates 
a binary random number, Ij g, which is one with probability pj, 
i.e., p[Ij^i = 1 ) = 1 - p{Ij^i = 0 ) = pj. If this random number is 
one, then Uj generates an ALC coded symbol, and sends 
it to the BS. The BS also generates the same binary random 
number as the one in Uj as they have shared the same random 
seed, thus it already knows in which time instants device Uj 
is sending a coded symbol. Using this random seed, both Uj 
and the BS can construct the same ALC code structure, so 
the BS is able to perform the decoding on the received coded 
symbols. 

Lor simplicity, we assume that devices use the same code 
degree dc and the same weight set Ws to generate ALC coded 
symbols. Let U denote the set of devices which are currently 
active and transmitting ALC coded symbols to the BS based on 
their access probabilities in the same RB. Then, the received 
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Fig. 3. Original AFC code graph at the BS. 


signal at the BS in time instant ^ is given by 

yt=^ + 2 ^, (7) 

leu 


where is the AWGN with variance cr|. As can be seen in 
O, the received signal at the BS is the noisy version of the 
sum of coded symbols of devices in U. Since coded symbols 
of each device have been generated by a linear combination 
of modulated information symbols, the received signal at the 
BS can also be seen as a coded symbol of an analog fountain 
code, where modulated information symbols are from devices 
in U. By using ([3i, (|7]i can be rewritten as follows: 

k 

yi= hj Ij,l 

j&A r—1 


where is the modulated information symbol of Uj and 
is the respective weight coefficient. Fig. [3] and Fig. |4] show 
the original and equivalent bipartite graphs of the AFC codes at 
the BS, respectively, when there are only two devices in S and 
dc = 2 . As shown in these figures, the equivalent bipartite graph 
can be considered as a bipartite graph of an equivalent AFC 
code where the information symbols are from both users. More 
specifically in the first time instant where Ii i = / 2.1 = l, both 
devices transmit coded symbols to the BS. But, in the second 
and third time instant, only one of the devices transmits a 
coded symbols to the BS, i.e., /i 2 = / 2,3 = l /i 3 = / 2,2 = 0 
. This is because of the fact that the devices transmit coded 
symbols based on their access probabilities. In fact, a device 
Uj with access probability pj transmit a coded symbols with 
probability ; otherwise it remains silent. 

Let us define = h/gPbP. Then is a random 
variable with the following probability distribution: 


Pixp = 




P-idc 

2H7 


if s = 0 , 

if s e Ws or - s e Ws. 


Thus, the mean and variance of Xpl can be calculated as 
follows: 


= E[X^"^]= V ^(.-.) = o, 
<7-= = V ^,2^1 


= E[xfP]= Y. 


sGWs 


kD 




( 9 ) 

( 10 ) 


where Since in M2M communications the 

number of devices are very large and X^ ;’s are independent 
random variables, according to the central limit theorem 
ZljGW Ef=i Xjg will approximately follow a Gaussian distri¬ 


bution with a zero mean and the following variance: 

ay = dccr^ 'Y, Pj\^j\^ ■ ( 11 ) 

jeu 

In other words, the received signals at the BS are normally 
distributed around the origin with average power ay. As the 
received coded symbols at the BS can be represented as coded 
symbols of an equivalent AEG code, the standard BP decoding 
algorithm can be applied to jointly decode all active devices 
at the BS. 

IV. Analysis of the Rateless Multiple Access for 
M2M Communications 

In this section, we analyze the performance of the proposed 
MA-AEC scheme by using a density evolution based approach. 
These results will be used as a basis to optimize the proposed 
rateless multiple access for M2M communications. 

As the BS determines the allocation of RBs to the devices, 
we assume a general case, where the number of active devices 
allocated to each RB is modeled by a Poisson random variable 
with parameter A. We also assume that each active device 
belongs to a delay group with a uniform probability. This 
model enables us to have a fair comparison with the throughput 
limit obtained in El. Thus, in the rest of the paper, we only 
consider one RB and optimize the access probabilities for the 
proposed scheme to simultaneously satisfy delay requirements 
for all devices. . 


A. Optimum Coordinated Multiple Access 

Before analyzing the proposed multiple access scheme, let 
us first investigate the optimum coordinated multiple access, 
which will be used as a baseline for the proposed probabilistic 
multiple access AEC scheme. 

We assume that the base station has already identified the 
devices which are transmitting in the same RB. The BS 
also knows all the channel information and the devices are 
transmitting with the same power. As shown in lll3[l . the 
maximum common rate that can be achieved by N devices 
transmitting in a typical RB with ordered effective channel 
gains {hj}^ and reference SNR 7 is given by: 


Ttc = min - log 2 ( 1-F 7 ) bps/Hz, (12) 

jGZf j \ ) 

which is simply obtained from the MAC capacity region when 
the transmission rate of all the devices are the same. This can 
be used as an upper bound on the average system common rate 
for the random access and later we will show how the proposed 
MA-AEC code can approach this limit. More specifically. 
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when the devices perform power control in a way to have the 
same received SNR 7 q at the BS, the common rate is given 
by: 

Rc = ^ log 2 (1 + Af7o) bps/Hz. (13) 

Therefore, the average system throughput increases with the 
number of devices with a logarithmic slope. The maximum 
payload size Lopt which can be successfully transmitted by 
each device in a resource block of duration ts and bandwidth 
W is then given by: 

Lopt = log 2 (1 + Lf'ro) bits, (14) 

which is obtained without considering the overhead for the 
contention phase. In fact, (fl4l i provides an upper bound on 
the maximum payload size which can be delivered from each 
active device at the BS. 


B. Asymptotic Performance Analysis of Multiple Access AFC 
based on the BP Decoding 

A commonly used analytical tool for analyzing a BP de¬ 
coder is density evolution ll 28 n . which calculates the evolutions 
of the message passing in the iterative decoding process. In 
this paper, we focus on the density evolution analysis in the 
asymptotic case, when the number of variable and check nodes 
go to infinity. 

Let us refer to information symbols of device as Type-Xj 
variable nodes in the bipartite graph of MA-AFC codes, for 
i e Each received coded symbol at the BS may connect 
to modulated information symbols of various sets of devices. 
A simple way to represent coded symbols at the BS is to 
divide them into different types according to their connections 
to different sets of devices. For simplicity, we represent each 
non-empty set of devices by a vector v of dimension N, where 
its entry, vi, is 1 if belongs to the set; otherwise, it is 
zero. Here, we refer to coded symbols which are connected to 
the modulated information symbols of set v as Type-v check 
nodes. It is easy to show that the probability that a coded 
symbol is of Type-v, qy, is given by: 


N 

9v= n 

i=l 




( 15 ) 


and there are in total 2^-1 various types of check nodes 
due to the fact that the total number of nonempty subsets of 
devices is 2^ - 1. 

Since there are 2^ - 1 types of variable and check nodes, 
we need to analyze the message between each specific types 
of variable and check nodes. In an AND-OR tree 

analytical method was proposed to analyze the decoding error 
probability of LT codes 113 ill with different types of variable 
and check nodes for erasure channels. However, the AND-OR 
tree approach cannot be applied directly to wireless channels. 
Here, we extend the conventional density evolution analysis 
to the case that there are multiple types of variable and check 
nodes. 

Lemma 1: Let us consider an M2M system consisting of 
N MTC devices wanting to simultaneously transmit their 
messages to a common BS. Each device is assigned with 


an access probability pj and code degree dj. Let denote 

the message passed from a Type-Xj variable node to a Type-v 
check node in the iteration of the BP decoding algorithm 
at the BS. Then approximated by a normal 

distribution with mean and variance 2rnP\ where 
can be calculated as follows: 

=fhi\^<rwdvi ^ ^ 2 

Y J- H“ O’y 

Vi = l 

where = mdjk, al = ^ 

N 

0-2 = > (17) 

2^1 

y 2 

and n{x) = (1 - tanh(a: - y^/x)) e T dy. 

The proof of this lemma is provided in Appendix lAl The BER 
for the device after t iterations of the BP decoder, denoted 
by can then be calculated as follows: 



where Q{x) = e I'^dz. 

C. Discussions on the asymptotic MA-AFC decoding perfor¬ 
mance at high SNR 

As stated before, in MA-AFC, the BS receives the sum 
of the coded symbols generated from active users. The BS 
thus tries to decode each device’s data from the received 
signals by performing the BP decoding. Now let us consider 
the asymptotic performance of MA-AFC as SNR goes to 
infinity, thus ignoring the noise. The decoding problem in 
this case is equivalent to solving a system of binary linear 
equations, where the variables are the information symbols of 
all active devices. More specifically, the BS requires to solve 
the following system of binary linear equations to find fepl for 
j e and l e zj, given the value of yi and matrices 

N k 

= ( 19 ) 

j=ll=l 

Let us define matrix G = [G^^1|G^^1|...|G^^1] and binary 
vector b = It is clear that G is the generator 

matrix of the equivalent AFC code at the BS and we have 
Y = Gb^ Since we consider that different devices select the 
weight coefficients from the same set Ws, it is possible that two 
different devices select the same weight coefficients for their 
transmitted symbols at the same time. More specifically, at 
time instant i, the probability that coded symbols of device 
and Um have at least one common weight coefficient, denoted 
by Pi^m, is given by 



where D is the weight set size and we assume that D d£ and 
D > dm. It is clear that when at least two columns of G are 
exactly the same, then ( fT^ does not have a unique solution. 
The following lemma gives the probability that at least two 
column of matrix G are exactly the same. 
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overhead can be calculated as follows; 



Lov 


— Lcp,l + ^cp,2 + ^cp,3 

= Lcp,l + Ns riog2 {N/Ns)^ 


+ Ns 


log2 



1 _ (1 _ e)(JV/iVs)-l 


( 22 ) 


As the BS allocates RBs to devices and devices with the 
same RA preamble will transmit at the same RB, each device 
then knows at which RB it should transmit its message based 
on its delay requirement. In this case, the amount of overhead 
can be further reduced to 


Fig. 5. The probability that at least two columns of the generator matrix at 
BS are the same versus the ratio of the number of coded symbols and that of 
information symbols for different weight set sizes. Each device has the access 
probability of 0.1 and code degree of 4. 


Lemma 2: Let us consider that m coded symbols are re¬ 
ceived at the BS and each information symbol of each device 
is connected to at least one coded symbol at the destination, 
i.e., the set of non-zero elements in each column of G is not 
empty. Moreover, we assume that different devices have the 
same access probability p and code degree dc. Let q denote 
the probability that two columns of G are exactly the same, 
then q can be calculated as follows: 


9 = 


(1+/3— dv) {tti —dij + l) I 
- ^ 


J^dqj — 1 




( 21 ) 


where jS = mpddk and dv = T/?]. 

The proof of this lemma is provided in Appendix |B] From 
this lemma, as can be also seen in Fig. |5l even when the 
weight set size is small, we can always uniquely decode 
all devices’ information symbols by transmitting more coded 
symbols. Moreover, when D is small, devices are required 
to transmit more coded symbols to guarantee that the BS 
can uniquely decode all devices information symbols, which 
clearly decrease the overall system throughput. To achieve 
higher throughput one may consider larger code degree dc, 
which on the other side significantly increases the decoding 
complexity at the BS that exponentially increases with dc. 


D. How much overhead is required in the contention phase? 

In the first phase of the contention phase, each active 
device sends a RA preamble to the BS. Let Lcp,i denote 
the length of the RA preamble. In the second phase, the BS 
sends the information about the number of devices selected a 
particular RA preamble. The average number of devices which 
have selected the same RA preamble is N/Ns, where N is 
the number of active devices and Ns is the number of RA 
preambles. These information can be broadcasted by the BS 
by using at least Lcp ,2 = A^sriog 2 (N/Ns)] information bits. In 
the third phase of the contention phase, the users will select 
a random seed of length L(s) which is calculated by (|6]l, for 
a given value of e and a given RA preamble s. Thus, at most 
Lj,p ,3 = NsL{s) information bits have to be sent to the BS in 
the third phase of the contention period. Therefore, the total 


Lov — Lcp^i -h riog2 (N/Ns)] 


log2 1 + 


1 _ (1 _ e) (N/Ns)-1 


(23) 


Let k denote the message length of each device and Tic 
denote the achievable common rate for a given arrival rate 
A using the proposed probabilistic multiple access approach. 
Then, this message can be successfully transmitted to the BS 
within a resource block of duration ts and bandwidth W if the 
following condition is satisfied; 


k Lov ^ tsM^T^c. 


(24) 


In other words, the maximum message size (in bits) which 
can be supported by the proposed approach in an RB can be 
approximated by; 

kmax = max{ rsVFRc — Lov,0}. (25) 


V. Optimization of Access Probabilities With 
DELAY Guarantees 


In this section, we formulate an optimization problem to si¬ 
multaneously satisfy the delay requirements of all the devices. 
The access probability and code degree for each device will 
be optimized in a way that delay constraints of all devices 
are satisfied. Here we minimize the average degree of the 
received signal at the BS to further reduce the complexity of 
the decoder. Therefore, the BS should perform the following 
optimization process: 

N 

min ^ pdi 

i=l 


subject to (a). Q < 5, (b). l < di < D, and (c). 

0 < Pi < for * s where 5 > 0 is the target bit error rate, 
m(°°) = lim rrS^\ and 

^ t—>-oo ^ 


it) Tl^iPo-^ciiL V- 

m- = - > 

i JU 




(t-i) 
a . 

3 


qy. 


It is important to note that Q.{x) is very close to zero for 
large values of x, so we can approximate ( fThl l as follows: 


(t) 27|hipcr2 

m - 5^ - 

^ k 


(26) 


for large x values. This means that the device Um which has 
the highest value of \hm\'^dmpm, will be decoded sooner than 
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Fig. 6. BER of device U 2 versus the BER of device Ui for the case that 
both the users have the complete access to the channel and di = d 2 . 


other devices at the BS. Let Pg ^ denote the bit error rate of 
device Ui at the BS. Thus, by using (l26T l we have; 


J 


\hj\‘^djpj ’ 


(27) 


where we assume that is constant which is valid when the 
number of devices is very large. Therefore, according to (fTSl l 
we have; 


pit) 


e,i 



\hi\'^diPi \ 
\h^\^d^Pj ) 


(28) 


More specifically, when the devices perform power control to 
have the same received SNR at the BS and use the same code 
degree dc, then (l28T l can be further simplihed as follows; 



which clearly shows that the device which has a higher access 
probability is decoded sooner than those with lower access 
probabilities. Fig. | 6 ] shows the BER for device U 2 versus the 
BER of device Ui for the case that both users have complete 
access to the channel, i.e., p^ = p 2 = 1 . As can be seen in 
this hgure, when the channel gains are the same, the BER 
performance is the same. Moreover, when the ratio of the 
channel gain of Ui and that of U 2 decreases, the BER of 
U 2 also decreases. 

By using ( |29] ), we can simplify the optimization problem 
as follows, assuming that all the devices use the same code 
degree dc and perform the power control to have the same 
received SNR 70 at the BS. 


N 

min ^ Pi 
i=l 


subject to (a), pi > and (b). 0 < < 1 , for i e zf. 

Therefore, the optimized access probabilities can be found as 
follows; 


Pi,opt 


min 


‘idc'jQti 


(30) 


According to ( |30] | we have Pi^opt l/L- Thus, if device 
has the access probability of pj, then the access probability 



Fig. 7. Physical layer configuration in LTE-A l23l . 

for Uj will be simply pj = ti/tj, which further simplifies the 
optimization process at the BS. Eurthermore, as the devices 
which have selected the same RA preamble will transmit 
in the same RB with the same access probability, the BS 
does not need to calculate the optimal access probability for 
every device. Instead, it calculated the access probabilities for 
each group of devices which have selected the same RB, and 
allocate the same access probability to devices within the same 
delay group. 


VI. Simulation Results 


In this section, we investigate the performance of the 
proposed MA-AEC scheme in M2M communication systems. 
We hrst compare our proposed MA-AEC approach with the 
conventional massive RA technique in the LTE-A standard 
a. Like lEill . as shown in Eig. [T] we assume that time is 
divided into time frames of length 10 msec and the total system 
bandwidth is 20 MHz. Each frame is further divided into 10 
subframes of length 1 msec. A subframe consists of 100 RBs, 
each with bandwidth W = 0.2 MHz. We also assume a total 
number of Ns = 60 preambles and for simplicity we assume 
several configurations of RAs, which includes 1 RA attempt 
per two time frames, 1 RA access attempt per time frame, 
and 2 RA attempts per time frame. Each successful device 
will then be allocated with 16 RBs for the data transmission 
in next subframes. We assume that the packet length of devices 
is small (few bytes), thus a device can successfully transmit its 
packet within one time frame. The conventional random access 
approach, which is currently used in the LTE-A standard uses 
backoff approach H. In this approach, when a random access 
attempt is unsuccessful, the device will wait for a random time 
before transmitting the new access request. More specifically, 
before the random access request, the device will wait for 
time where is randomly drawn from range [0,Xi - 1 ]. 
The device can send a maximum number of Nra RA request 
using this backoff approach. More specihcally, we assume 
maximum number of RA attempts Nra = 10 and backoff in¬ 
dexes Xi e {10, 20,30,40,60,80,120,160, 240,320,480,960} iS . 

Eor the comparison, we use two metrics, blocking proba¬ 
bility and average access delay. The blocking probability is 
dehned as the ratio of the MTC devices failing to access the 
BS due to exceeding their maximum allowed random access 
attempts to the total number of MTC devices. The average 
access delay is defined as the average time elapsed from the 
time instant when a device sends the first access request until 
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Fig. 8. The blocking probability for MTC devices versus the number of 
devices. The total number of preambles is Ns = 60. All device in the MA- 
AFC case use the same code degree 8. 



Fig. 9. The average access delay for MTC devices versus the number of 
devices. The total number of preambles is Ns = 60. The access probability 
for MTC devices in Scenario 1 is Ns/N and in Scenario 2 is lONs/N. All 
device in the MA-AFC case use the same code degree 8. 

that it succeeds. Fig. 0 shows the blocking probability for both 
the proposed and conventional approaches. As can be seen 
in this hgure, the blocking probability for the conventional 
approach is very high especially when the number of devices 
is large. However, devices in the proposed scheme can access 
to the BS with a very small blocking probability, due to the 
fact that the devices can transmit simultaneously and they do 
not need to wait until they select separate preambles. 

Fig. |9] shows the average access delay versus the number of 
devices. As can be seen in this hgure, the proposed approach 
outperforms the conventional approach and have roughly 70% 
less access delay. In Fig. |9] we considered two scenarios for 
MA-AFC. The access probability for all devices in Scenario 
1 and Scenario 2 is set to Ns/N and lONs/N, respectively. It 
is important to note that the proposed approach in Scenario 2 
has a lower access delay due to the fact that more devices can 
simultaneously transmit but it has higher complexity at the BS 
due to larger degrees of equivalent AFC code at the BS. 

We then simulate our proposed approach within one RB, and 
compare the results with the fundamental throughput limits in 
10 . For this aim, we assume that the reference SNR is 7 = 0 
dB. This is equivalent to the case that a device transmits with 
10 dBm (lOmW) power over IMHz bandwidth, a noise power 


Fig. 10. Maximum achievable common rate versus the arrival rate for the 
case that the devices perform power control to have the same received SNR 
of 0 dB at the BS. 

spectral density of -174dBm/Hz, a receiver noise hgure of 5dB, 
a receiver antenna gain of 14dB, a 3.76 pathloss exponent, a 
128dB pathloss intercept at 1000m, and a cell radius of 1360m 
1 131]. We also assume that ts = 1 sec. 

Let us hrst investigate the maximum common rate which 
can be achieved by the proposed probabilistic MA-AFC ap¬ 
proach in the contention-free case. Fig. [TO] shows the maxi¬ 
mum common rate versus the arrival rate for the case that the 
devices perform the power control to have the same received 
SNR of 0 dB at the BS. We consider three cases for the 
MA-AFC approach with different access probabilities. As can 
be seen in this hgure, the achievable rate of the proposed 
approach is very close to the optimal coordinated multiple 
access 0. It is important to note that in the proposed 
approach we design the access probabilities in a way that 
the number of devices which simultaneously transmits is 4, 
6 , and 8, respectively, in order to minimize the decoding 
complexity at the BS. Fig. [TO] shows that even with a small 
access probability, the achievable common rate of the proposed 
approach is very close to the optimal coordinated approach, 
especially when the arrival rate is very large. 

In Fig. [TT] we show the maximum payload size which 
can be supported by the proposed approach in a resource 
block of duration ts = 1 sec and bandwidth w = 10 kHz 
for different arrival rates. An upper bound on the maximum 
payload size which has been obtained from (fl4li is also shown 
in this hgure. As can be seen in this hgure, as the arrival rate 
increases, the payload size which can be supported decreases. 
This is because that the common rate is a decreasing function 
of the arrival rate (Fig. [loll. Moreover, as the arrival rate 
increases, the overhead calculated in (|2^ increases, which 
reduces the supported payload size. Also, by increasing the 
access probability, a higher common rate can be achieved, 
which leads to higher supported payload sizes. 

Fig. [ 12 ] shows the average BER versus the number of 
received coded symbols at the BS when the code degree for all 
devices is dc = 8 and the reference SNR is 7 = 30 dB. As can be 
seen in this figure, by increasing the number of devices, the BS 
requires more coded symbols, in order to decode all devices’ 
information symbols. Moreover, by increasing the number 
of devices, the average number of required coded symbols 
per device remains relatively constant to achieve a certain 
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Fig. 11. Maximum supported payload size versus the arrival rate. Devices 
perform power control to have the same received SNR of 0 dB at the BS. 
The devices use the same code degree of 8 and the same weight set Wg. 



Fig. 12. BER versus the ratio of the number of coded symbols and the total 
number of information symbols, m is the number of coded symbols, N is 
the number of devices, and k is the length of the message for each device. 

BER. This indicates that the proposed scheme can provide 
a certain level of BER at the BS for each device even with a 
large number of devices and the decoding performance is not 
significantly degraded by increasing the number of devices. 

We now consider the case that devices have different delay 
requirements. We assume that the payload size of each device 
is 1 kb and 100 MTC devices belongs to each specific 
delay group, as shown in Table [III By using the optimization 
approach in Section IV, we can derive the optimized access 
probability for each group. Eig. [T3 shows the average BER 
versus the delay for devices in each group of MTC devices 
specified in Table [III when all the devices has the same 
code degree dc = 8. As can be seen in this figure, the 
delay requirement for devices in various groups have been 
simultaneously satisfied. 

VII. Conclusions 

In this paper, we proposed a novel rateless multiple access 
scheme to support a large number of devices in M2M com¬ 
munications. This is achieved by representing the multiple 
access process by a capacity approaching analog fountain 
code to minimize the access delay while achieving a higher 
throughput. In the proposed scheme, devices are allowed 


Fig. 13. The average bit error rate for different groups of devices classified 
by delay requirements. The total number of devices is = 1000, each 
group has 100 devices; the devices use the same code degree dc = 8, access 
probabilities are calculated from (30) and 7 = 0 dB. Dashed lines show the 
delay requirement for different groups of devices. 

to transmit at the same time and the same channel, thus 
forming an equivalent analog fountain code at the base station 
(BS). The BS then performs joint decoding to recover all 
devices’ messages. We further analyzed the proposed multiple 
access analog fountain code by using the density evolution 
technique and formulated an optimization problem to find the 
optimum access probability for each device to satisfy the delay 
requirements of MTC devices. Simulation results showed that 
the proposed scheme closely approaches the optimal common- 
rate performance of the multiple access channel, even for a 
very large number of devices. Moreover, the proposed scheme 
can satisfy the delay requirements for all the devices with a 
relatively low decoding complexity, which increases linearly 
with the number of devices. We also proposed a simple 
contention-based strategy and formulated the average overhead 
of the proposed scheme. Our novel random access and rateless 
multiple access schemes also achieved a significantly lower 
access delay and blocking probability compared to existing 
massive random access schemes for cellular-based M2M com¬ 
munications. 


Appendix A 
Proof of Lemma[T| 


Let us consider the equivalent AEC code at the BS, where 
the equivalent generator matrix and the message vector are 
denoted by G and b, respectively. We first define Yj = 
Y.r'^M(j)\r93,r'^r'^ where M{j)\r is the set of all variable 
nodes connected to check node yj except variable node hr. In 
S, it has been shown that the average log likelihood ratio 
(LLR) passed from check node yj to variable node br in the 
iteration of the BP algorithm, can be approximated 

as follows: 


:|lW ] = 2g,r r°°^dy 


/— cso fo{y) 

where fm (y) is defined as follows for m e {0,1}: 


fm{y) = XI 

r'^M{j)\r 


y/27r 


{y-Yj)"^ ]7 PxiK')- 

r'eM{j)\r 
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As shown in when the number of variable and check 
nodes go to inhnity, /o (y) i ^ ^ 

is the variance of Yj. As shown in S, in the asymptotic 
case the message passed from a variable to a check node is 
normally distributed and we can hnd the mean and variance of 
the message passed from a Type-Xj variable node to a Type-v 
check node, denoted by , , as follows; 


+ (7^ 


E 


ie{4U}i = t 

^ e'eJ^{r)\e 

var{4;^J=2|E{LW Jl, 




where is the set of all check nodes connected to a 

Type-Xi variable node except check node As we assume 
that the number of variable and check nodes go to infinity, 
according to the law of large number, for large code degree d^, 
J2g'^U{Xi)\£9"ii X approximated by a'^mdi/k, where 

m is the number of coded symbols. As a Type-Xj variable 
node is connected to a Type-v check node with probability gv, 
then the average total LLR value for a Type-Xj can be found 
as follow; 


='1\hi\'^a■l,dv^ ^ ^ 3 gv, 

\T 1- f7v 


Vi = l 


where o-y = ) and f2(a;) = 

4= (l - tanh(x - y^/x)) e~^dy This completes the 

proof. 


Appendix B 
Proof of Lemma|2] 

To generate each AFC coded symbol, dc information sym¬ 
bols are selected from those with minimum degrees, so after 
receiving m coded symbols at the BS, the degree of each infor¬ 
mation symbol will be either dv and dv - 1, with probabilities 
l+p-dv and dv-0, respectively. Moreover, two columns of G 
are the same if the number of their nonzero elements are the 
same. The probabilities that two columns of G have the same 
degree dv and dv-l, are (l+p-dv)"^ and {dv-0)‘^, respectively. 
Furthermore, the probability that the nonzero elements of the 
two columns of the same degree are the same is l/f jD“, 

where a is either dy or dy - 1. This arises from the fact that 
each nonzero element of each column is randomly chosen from 
m possible positions, and the value of each nonzero element 
is randomly selected from D possible weights. This completes 
the proof. 
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